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I .  INTRODUCTION 

There  have  been  two  kinds  of  methods  published  for  creating  the 
ship/wave  problem  for  a  ship  moving  at  very  low  speed: 

(1)  First  Ogilvie  (1968)  used  an  order-of-magnitude  argument  to 

obtain  a  linear  free-surface  condition  that  would  lead  to  the  prediction 
of  a  plausible  wave  motion  at  very  low  speeds.  There  were  two  essential 
points:  (a)  the  waves  should  have  very  short  wavelength,  and  (b)  the 

waves  should  propagate  on  the  nonlinear  streaming  flow  around  the  correspond¬ 
ing  double  body.  In  the  linear  wave  problem,  the  apparent  cause  of  wave 
generation  is  an  effective  pressure  distribution  on  the  free  surface,  which 
arises  mathematically  because  the  double-body  flow  does  not  really  satisfy 
the  precise  free-surface  conditions.  Ogilvie  treated  only  the  case  of  a 
fully  submerged  two-dimensional  body.  Later,  Baba  and  Takekuma  (1975)  extended 
this  concept  to  solve  the  problem  of  a  three-dimensional  surface-piercing  body, 
e.g. ,  a  ship.  They  went  so  far  as  to  derive  a  wave-resistance  formula  based 

on  this  approach.  Maruo  and  Fukazawa  (1979)  extended  this  approach  further, 
using  a  coordinate  transformation  to  simplify  the  analysis. 

(2)  Keller  (1974)  developed  the  first  ray  theory  for  the  low-speed 
problem.  Inui  and  Kajitani  (1977)  used  a  procedure  based  on  a  method  of 
Ursell's  (1960),  which  is  essentially  a  ray  method.  Later,  Keller  (1979) 
further  developed  his  ray  method  with  systematic  asymptotic  expansions,  and 
he  applied  his  theory  to  a  thin  ship  and  a  special  class  of  "streamlined" 
ships.  In  Keller's  ray  theory,  the  waves  are  apparently  generated  only  at 
the  stagnation  points  on  the  body;  the  amplitude  and  phase  of  the  waves 
are  then  modified  gradually  by  the  nonuniform  flow  around  the  double  body. 

There  are  difficulties  in  both  methods.  In  the  first,  the  linear 
free-surface  condition  can  be  written  in  such  a  way  that  the  terms  on,  say, 
the  left-hand  side  are  rapidly  varying  wavelike  quantities,  while  those  on 
the  right-hand  side  are  slowly  varying  in  space.  The  latter  are  completely 
known;  they  represent  the  fictitious  pressure  field  imposed  on  the  free 
surface.  The  situation  can  be  compared  to  the  much  simpler  problem  cited 
by  Keller  (1979)  (see  his  Appendix  A) : 
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u" (x)  +  k‘u (x)  =  q (x)  , 

where  g(x)  is  analogous  to  the  fictitious  pressure  distribution.  If  k 
is  very  large,  wave  solutions  of  this  differential  equation  reprer.-  :.L  very 
short  waves.  The  general  solution  of  the  above  equation  can,  of  course ,  be 
written  out  explicitly,  and  it  can  then  be  expanded  asymptotically  ft  • 
k  -*•  00  .  If  the  domain  of  x  is  <  x  v  +®  ,  the-  only  part  of  the  as/nu  t.  ...tic 

expansion  that  represents  waves  comes  fr  in  the  homogeneous  problem,  -.nd  1  t.s 
amplitude  and  phase  can  be  determined  g  r.i  ral  ly  only  if  som>hcw  they  are- 
known  at  some  point,  possibly  at  infinity.  In  addition,  there  is  a  part icul ab¬ 
solution  ,  which  can  be  represented  asymptotically  as  a  series  in  inverse 
powers  of  k  ;  it  represents  a  slowly  varying  solution  if  g(x)  is 
slowly  varying  (as  assumed).  If  the  domain  is  restricted  to,  say,  0  •  x  <"  , 
the  wave  part  of  the  solution  depends  entirely  on  the  values  of  q(x)  and 
its  derivatives  at  x  =  0  and  on  the  two  boundary  conditions  imposed  on 
u(x)  .  So  we  can  say  that  the  generation  of  waves  in  such  a  case  is  unaffected 
by  the  function  g(x)  except  in  a  neighborhood  of  x  =  0.  This  raises 
a  doubt  about  the  fundamental  supposition  of  the  first  method,  namely,  that 
the  waves  are  generated  (mathematically  speaking)  by  the  fictitious  pressure 
distribution  on  the  free  surface.  All  that  matters  is  the  behavior  of  g(x) 
near  x  =  0. 

Dagan  (1972)  pointed  out  that  the  base  flow  on  which  the  waves  propagate 
should  not  be  just  the  double-body  flow,  as  assumed  by  Ogilvie  and  others,  but 
at  least  two  terms  in  the  "naive  expansion,"  which  is  the  expansion  that  is 
obtained  if  the  problem  is  expanded  formally  and  strictly  in  terms  of  power 
series  in  the  Froude  number.  Keller  (1979)  arrived  at  the  same  conclusion 
and  noted  further  that  then  the  fictitious  pressure  distribution  vanishes  from 
the  wave  problem.  That  is,  the  free-surface  condition  becomes  homogeneous. 

In  fact,  in  Keller's  (1979)  ray  theory,  the  wave  part  of  the  velocity 
potential  function  satisfies  the  Laplace  equation,  which  is  homogeneous,  as 
well  as  homogeneous  free-surface  and  body  boundary  conditions.  Thus  any 
solution  that  is  found  can  be  multiplied  by  an  arbitrary  constant.  Keller 
introduced  a  so-called  "excitation  coefficient"  for  certain  simple  special 
cases.  Still,  his  method  fails  near  the  stagnation  point  of  the  double-body 
flow:  From  the  dispersion  relation,  the  wave  number  becomes  inifinte  there. 


and  the  amplitude  of  the  wav-  p  becomes  infinit  coo.  Since,  as  Kel  i-.r  hirnseif 
pointed  out,  the  generation  of  waves  in  the  short-wav* •  problem  depends  essentially 
on  conditions  at  the  boundary  point  (consider  again  the  simple  differential- 
equation  problem  cited  above),  the  failure  of  tin.-  ray- theory  assumptions  near 
a  stagnation  point  seems  to  be  crucial. 

We  present  here  the  second  part  of  a  study  to  resolve  these  quest ions . 

The  f' rst  part  is  reported  by  Ogilvie  and  Chen  (1982);  it  will  be  referred  to 
subsequently  simply  as  "I".  They  d*  rived  .*  nonhomogeneous  body  boundary  condition 
for  the  wave  part  of  the  potential  function.  Their  free-surface  condition  is 
homogeneous,  as  required  from  the  work  of  Keller.  The  nonhomogeneity  of  the 
body  boundary  condition  is,  as  we  shall  show,  adequate  for  determinina  the 
solution  without  any  arbitrary  additive  or  multiplicative  quantities. 

We  use  the  method  of  matched  asymptotic  expansions  to  solve  the  low-speed 
problem  in  two  dimensions.  Our  outer  region  is  a  thin  layer  near  the  free 
surface,  far  behind  the  body  (in  terms  of  wavelength) .  We  use  the  generalized 
WKB  method  to  determine  the  nature  of  the  wave  motion  in  this  region;  this 
method  is  very  similar  to  a  ray  method.  Then  we  formulate  a  near-field 
problem,  applicable  to  a  very  small  region  near  the  stagnation  point.  Con¬ 
siderations  of  the  flow  properties  near  a  stagnation  point  are  found  to  be 
sufficient  for  determining  the  near-field  solution  completely.  Then  matching 
to  the  outer  solution  permits  the  latter  to  be  determined  as  well. 

The  final  formula  for  the  wave  resistance.  Equation  (89),  is  surprising: 

Wave  resistance  is  proportional  to  »  where  U  is  the  forward  speed.  We 

believe  that  this  is  the  correct  asymptotic  relationship,  although  it  is  not 
likely  to  be  useful  to  a  naval  architect.  It  clearly  has  no  range  of  validity 
in  U  in  which  its  predictions  overlap  those  from  conventional  wave-resistance 
analyses.  What  is  still  needed  is  a  small-U  solution  that  gives  this  formula 
as  U  -*•  0  and  gives  the  results  of  conventional  linear  thoery  as  U  -+  00  . 

There  is  an  even  stronger  contrast  between  the  low-speed  theory  and 
conventional  linear  theory  in  the  case  of  a  submerged  body.  (The  low-speed 
theory  predicts  that  wave  resistance  is  proportional  to  e  as  U  -*  0  .  ) 


Tulin  (1982)  has  produced  an  exact  theory  that  bridges  f  twc 
much  more  difficult  to  find  a  comparably  general  theory  for  tl 
surface-piercing  body. 


II.  OUTER  SOLUTION  BY  THE  »K3  METHOD 


We  want  to  find  a  velocity  potential  LU?(x,y)  for  the  streaming  flow  past 
a  two-dimensional  body  that  intersects  the  free  surface.  The  stream  has  speed 
rJ  in  the  positive-x  direction.  The  undisturbed  tree  surface  lies  on  the  x 
axis,  with  the  y  axis  directed  upwards.  We  take  uhe  origin  of  coordinatr-s 
inside  the  body  in  such  a  way  that  the  intersections  of  the  body  and  the 
undisturbed  free  surface  are  located  at  '-xOl0)  and  (xg,0)  .  The  latter  is 

the  one  of  primary  interest  in  this  paper,  since  it  lies  on  the  downstream  side 
of  the  body.  All  length  dimensions  have  been  normalized  with  respect  to  L  , 
any  convenient  characteristic  length  of  the  body.  The  small  parameter  of  tie 
problem  is  taken  as 

£  =  F2  =  U2/gL  ,  (1) 

where  g  is  the  gravitational  constant  and  F  is  a  Froude  number.  The  shape 
of  the  free  surface  is  given  by  a  relationship  y  =  eH(x)  ,  where  H(x)  is  to 
be  determined  as  part  of  the  solution  of  the  problem. 

The  statement  of  the  problem  is  as  follows: 
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(7) 

These  are  the  same  as  in  (I) ,  although  the  [K]  condition  was  not  explicitly 
used  there.  Either  the  [F]  condition  or  the  [K]  condition  is  redundant. 

We  assume  that  the  solution  can  be  divided  into  two  parts : 


i'[>  ,y) 


(8) 


“  - 

=  +  K{x,y;r.,' 

H  (x)  =  Mix;,' )  +  Mix;,)  .  ( ■)) 

The  first  part,  represented  by  v  ( x ,  y ;  r. )  and  H(x;<.)  ,  is  the  so-called 

expansion  (see  (I));  it  is  the  formal  solution  that  is  obtained  by  simply  sub¬ 
stituting  a  power  series  in  e  into  the-  conditions  (2)  -  (7).  It  does  n  ,t 
represent  a  wavelike  motion,  although  there  is  a  corresponding  free-surface 
deformation  near  the  body,  which  vanishe-  downstream  as  well  as  upstream.  The 
other  terms  in  (8)  and  (9)  represent  true  wave  motions  (suggested  by  the  nota¬ 
tion  ~  ) .  We  further  assume  that  the  wave  part  of  the  solution  can  be 
expanded  as  follows: 

4»(x,y;r.)  't  eot+1<j>1  (x,y;  o,Y)  +  ea+2(J2  (x,y;  Q,Y)  +  ...  ; 

H (x; c)  ^  Canj(x;0)  +  ta+1n2(x;0)  +  ...  . 

Three  new  quantities  have  been  introduced: 

(i)  a  is  a  real  number  greater  than  unity.  In  (I) ,  it  was  taken  as 
1  .  Actually,  that  is  simply  the  smallest  possible  value  of  a 
that  leads  to  a  linear  problem  for  the  wave  motion.  Now  we  must 
use  the  nonhomogeneous  boundary  condition  developed  in  (I)  to 
determine  the  correct  value  of  a  .  Note  that  the  difference  in 
the  powers  of  e  between  (10)  and  (11)  results  from  the  [H]  con¬ 
dition,  (3)  . 

(ii)  0(x)  is  a  rapidly  varying  phase  function,  which  we  shall  also 
write  in  the  form 

0 (x)  =  0(x)/e  .  (12) 

We  shall  assume  that  0(x)  is  slowly  varying  in  the  sense  that 
its  derivative  is  of  the  same  order  of  magnitude  as  0  itself. 

(iii)  Y  is  a  stretched  coordinate: 

Y  =  y/e  .  (13) 

Effectively,  we  treat  this  as  a  multiple-scale  problem,  x  and  y  being  used 
to  describe  changes  that  occur  on  a  scale  comparable  to  the  dimensions  of  the 
body,  0  and  Y  being  used  to  describe  the  details  on  the  scale  of  the  wave¬ 
length.  We  imply  here  that  the  generated  waves  have  wavelength  that  is  0(e)  . 
This  is  valid  in  the  outer  region,  which  means  a  region  many  wavelengths  away 
from  the  downstream  stagnation  point.  The  implication  is  not  valid  very  near 
to  the  stagnation  point.  The  latter  case  is  discussed  in  the  next  section. 
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The  expansions  (10)  and  (11)  are  called  ger  aalised  WKB  expansions .  The 
functions  <p±  and  rj all  represent  wav  motions  that  are  superposed  on  a 
nonuniform,  nonwavelike  base  flow  given  by  $  and  H  . 


Let  us  define 


K(x) 


d9(x) 

dx 


We  note  the 

9$ 

9x 

3£ 

9y 

dH 

dx 


following  formulas  for  dif  f erer.tiations  : 

—  +  ea+  {(j> lx  +  [K/e] $1q}  +  ea+  {i2x  +  lK/e] <f>20l  +  •••  > 

f“  +  £a+1{^iy  +  [Ve]^ly}  +  ca+‘M^2y  +  tVe]i2x}  +  ...  ; 

+  £“{nlx  +  tK/e]nl0}  +  £a+1{n2x  +  tK/e)  n2Q}  +  ...  . 


(14) 


(15) 

(16) 
(17) 


Now  we  substitute  (8) -(9)  and  (10) -(11)  into  the  conditions  of  the  problem, 
starting  with  the  Laplace  equation  (2).  Then  we  rearrange  terms  according  to 
powers  of  e  and  set  the  coefficient  of  each  power  separately  equal  to  zero. 
From  the  coefficient  of  ea-1  ,  we  obtain: 

[L]  K2^l0Q  +  0lyY  =  0  .  (18) 

Similarly,  from  the  coefficient  of  ea  ,  we  obtain: 

M  k2*200  +  *2YY  =  -  {(K*10>X  +  K*10X  +  2$lyY}  ’  (19) 

The  letter  subscripts  indicate  partial  or  total  derivatives,  as  appropriate. 

(For  example,  Kx  =  dK/dx  ,  whereas  <j>x  =  9<}>/9x  .  ) 


Next  we  substitute  into  the  [H]  condition,  (3) .  Initially,  all  functions 
of  y  and  of  Y  must  be  evaluated  on  eH(x)  and  on  H(x)  ,  respectively. 

Then  we  expand  these  functions  in  Taylor  series  as  follows:  $  (and  its  deri¬ 
vatives)  is  expanded  with  respect  to  y  =  eh  ,  whereas  (jj^  is  expanded  with 
respect  to  y  =  0  and  Y  =  H  .  (The  last  is  permissible  if,  as  assumed, 
a  >  1  .)  We  note  that,  from  the  definition  of  the  naive  expansion, 

(x,eH)  =  eHx$x(x,eH)  .  (20) 


From  the  coefficient  of  e'~  in  (3),  we  obtain: 


[H]  rijUjG)  =  -  K(x)  $x(x,cHHly(x,0;0,H)  .  (21) 

Similarly,  from  the  coefficient  of  ra+1  , 

[H]  n2  =  -  $x{K$20  +  $ix  +  K5^l0y  +  Mly}  on  y  =  0  ,  Y  =  H  .  (22} 

Here  and  throughout  this  section,  the  notation  "  y  =  0  " ,  as  in  (22),  refers 
only  to  ij>^  .  As  already  mentioned,  and  as  indicated  explicitly  in  (21)  ,  we 
evaluate  $  on  y  =  eH  . 

Following  the  same  procedure  with  the  [K]  condition,  (4),  we  obtain; 


[K] 

*1Y  = 

Kixnl0 

on  y  =  0  ,  Y  =  H  ; 

ry  =  0  , 

(23) 

[K] 

*ly  +  *2Y  +  = 

ix  + 

K«x<i>iG  +  K+x^e  "  ^iyyH  on 

.Y  =  H  . 

(24) 

We  can  eliminate  rij 

between 

(21)  and  (23)  (or,  alternatively,  start 

with 

the 

[F]  condition,  (5)).  ' 

Then  the 

<t>j  problem  is  given  by  the 

following : 

[L] 

K2»100  +  *iyY 

-  o  —  2~ 

=  0 

in  y  <  0  ,  Y  <  H  ; 

(18) 

[F] 

+  1Y  +  K  ***l00 

=  0 

on  y  =  0  ,  Y  =  H  . 

(25) 

Then  we  eliminate  n2  between  (22)  and  (24) ,  which,  together  with  (19) ,  gives 
the  <J> 2  problem: 


[L] 

k2+200  *  *2YY 

=  -  {<K*l0>x  +  4lQx  +  20lyy}  in 

y  <  0  ,  Y  <  H  ; 

(19) 

[F] 

*2Y  ♦  rtfree 

=  "  ♦ly  +  *xnlX  +  “Sc*!©  '  H^lyY  " 

K*x*lxo 

-  K2*x5*ioey  -  K»xVlY0  - 

y  =  0  , 
on 

Y  =  H  . 

(26) 

Now  let  us  solve  the  problem  stated  in  (18)  and  (25) .  The  independent 
variables  are  0  and  Y  ;  we  can  consider  x  and  y  as  if  they  were  para¬ 
meters.  By  a  minor  change  of  variables,  we  transform  the  first  equation  into 
the  Laplace  equation.  Let 


Then  (18)  and  (25)  become: 


t.V  V 


<“  ' 


*1xx  +  =  0 

in  y  <  0  , 

Y  <  0  ,• 

(28) 

-2~ 

d> ,  a  +  <fc  <j>.  =  0 

Tiy  x^ixx 

on  y  =  0  , 

Y  =  0  . 

(29) 

The  slight  change  from  Y  to  Y  in  (25)  has  enabled  us  to  formulate  a  problem 
with  a  free-surface  condition  given  on  the  X  axis. 

We  introduce  complex  variables  Z  -  3  f(Z)  : 

Z  =  X  +  iY  ,  (30) 

^  =  Re{£ (Z) }  .  (31) 

(Of  course,  f  also  depends  on  x  and  y  ,  but  we  continue  to  treat  the 
latter  as  parameters.)  The  Laplace  equation,  (28),  is  automatically  satisfied 
when  (j^  is  defined  as  in  (31).  The  free-surface  condition,  (29),  can  now  be 
rewritten : 

Re{ i^(x,eH)f"  +  if'}  =  0  on  y  =  0  ,  Y  =  0  .  (32) 

Since  we  consider  x  (and  thus  H(x)  )  as  being  fixed  in  (32),  this  condition 
is  valid  for  all  X  .  Thus  we  can  define  a  new  function 

W(Z)  =  (Z)  +if'(Z)  ,  (33) 

which  can  be  continued  analytically  into  the  upper  half  plane  as  follows : 

W(Z)  =  -  W(Z)  ,  (34) 

where  the  bars  here  denote  complex  conjugates.  Since  W(Z)  is  analytic  in 
the  entire  lower  half  of  the  Z  plane,  we  now  conclude  that  it  is  analytic  in 
the  entire  plane.  Thus  it  must  be  a  constant.  From  (29) ,  its  real  part 
vanishes,  and,  without  loss  of  generality,  we  set  its  imaginary  part  equal 
to  zero  also.  Then 


$xf"(Z)  +  if ' (Z)  =  0  in  the  Z  plane. 


(35) 


This  is  an  ordinary  differential  equation  for  f  (Z)  ,  which  is  easily  solved: 


f  (Z)  =  A(x,y)  exp  {-iZ/t^.}  ,  ( V ) 

where  A(x,y)  is  an  arbitrary  constant  with  respect  to  X  and  V  .  In 
general,  an  additive  constant,  say  C(x,y)  ,  can  be  added  to  (36) ,  but.  it  con¬ 
tributes  nothing  to  the  wave  solution,  and  so  we  set  it  equal  to  zero. 

We  require  that  f  be  a  periodic  function  of  0  .  We  have  not  yot  spur i- 
fied  >•  =  eO(x)  ,  and  so,  without  loss  of  generality,  we  can  require  that  the 
period  be  2tt  .  Substituting  Z  =  U/K  s  xY  into  (36)  and  taking  into  account 
that  K  =  0' (x)  =  e G' (x)  ,  this  requirement  is  equivalent  to  the  following: 

0  „  9+2  tt 

-  — rr; +  2n  =  -  — =-3 - 

e$xO*  cix(0->2T\)' 

It  then  follows  that 


0  (x)  =  --  dC/**(C,eH(C))  ,  (37) 

e  x 

x0 

where  once  again,  without  loss  of  generality,  we  have  set  a  constant  of  integra 
tion  equal  to  zero,  substituting  into  (36)  and  then  using  (31) ,  we  have  the 
solution : 


4>j(x,y;0,Y)  =  i?e{f(Z)}  =  Re{ A(x ,y)  exp  ( (Y  -  H) /$£)  exp  (i0 (x)  )  }  .  (38) 


From  (21)  we  obtain  the  wave-elevation  function: 


whe  re 


hjtx;©)  =  /?e{b(x)  exp  (iO(x))}  , 


b  (x)  =  iA(x,0)/<fx(x,eH)  . 


From  now  on,  for  convenience,  we  drop  the  notation  Re  ,  but  we  imply  that  it 
should  be  included  in  expressions  like  (38)  and  (39). 

The  situation  represented  by  the  solution  and  n1  above  is  familiar 

in  applications  of  the  WKB  method:  We  now  know  the  basic  form  of  the  wave 
solution,  but  we  do  not  know  its  amplitude  anywhere.  We  know  the  phase  func¬ 
tion,  0(x)  ,  but  the  (complex]  amplitude  A(x,0)  is  completely  unknown. 


In  order  to  obtain  more  information  about  A(x,0)  ,  we  must  consider  the 
second-order  problem,  given  by  (19)  and  (26) .  The  form  of  the  left-hand  sides 


of  (19)  and  (26)  is  identical  to  that  of  (18)  and  (25)  ,  which  gave  the  first-order 
problem.  So  we  can  expect  the  second-order  solution  to  represent  waves  like  those 
described  in  (38)  and  (39).  That  is,  will  contain  a  part  that  satisfies  the 

homogeneous  counterpart  of  (25) ,  and  this  part  will  depend  on  the  same  phase  func¬ 
tion  as  that  in  ^  .  In  addition,  we  note  that  the  right-hand  sices  of  (ly)  and 
(26)  are  linear  in  <{q  and  n ^  ,  and  so  the  solution  of  the  nonhor.iogei.eouR  piv.n- 
lem  will  also  involve  the  same  phase  function.  So  we  now  write: 

<t>2  =  A2(x,y)  exp  {  (Y-H) ,  ■  x}  exp  { it  (x)  }  ,  (41) 

where  0(x)  is  still  given  by  (37)  and  A2(x,y)  is  an  unknown  function. 

When  (41)  is  substituted  into  (19)  and  (26) ,  it  is  evident  that  the  equations 
can  be  satisfied  only  if  the  right-hand  sides  of  those  equations  are  separately 
equal  to  zero.  This  provides  the  further  conditions  needed  for  determining  the 
first-order  solution.  Substituting  (38)  and  (39)  into  the  right-hand  sides  of 
(19)  and  (26) ,  setting  them  equal  to  zero,  and  letting  y  =  H  ,  we  find  that: 

iAKjj  +  2iKAx  -  2ikAHx/$x  +  2Ay/$x  =  0  ,  (42) 

^x^x  ~  ‘fyyb  “  2iHxA/4>x  -  i^xx'A/'J’x  =  0  ,  1 43) 

where  we  have  used  the  fact  that  Hx  =  _^x^xx  +  0^c^  •  which  follows  readily 

from  (3)  and  (4) .  We  use  this  relationship  again  in  (42)  ,  noting  also  that 

-2 

K  =  -  1/4X  ,  to  obtain  the  following: 

Ay(x,0)  =  iAx (x, 0)  .  (44) 

Then  we  use  this  result  and  the  definition  (40)  in  (43)  to  obtain: 

Ax  +  A^xx/^x  =  0  *  (45) 

This  differential  equation  is  easily  solved:  A(x,0)  =  Ag/4^  ,  where  AQ  is 
strictly  a  constant.  Thus  we  have  found  the  form  of  A(x,0)  ,  and  so  we  have 
for  ijjj  : 

4> L  (x,0;0,H)  =  (A0/4x(x,eH)  }  ei0^X^  .  (46) 

The  const?  AQ  can  only  be  found  from  matching  this  solution  with  a  near¬ 
field  (ir^o  )  solution.  This  situation  arises  because  our  fundamental  equation 
is  an  ellip  ic  equation  (the  Laplace  equation) ,  and  the  typical  WKP  wave  approxi¬ 
mation  is  not  valid  in  a  region  near  the  body,  where  elliptic  behavior  domi¬ 
nates  the  wave  behavior. 


III.  INNER  SOLUTION  BY  THE  SOURCE-DISTRIBUTION  METHOD 

The  solution  obtained  in  the  preceding  section  represents  a  vivo  motion 
with  very  short  waves.  To  be  precise,  from  (37)  and  (38),  it  is  evident  that 
the  local  wave  number  is 

G' (x)  =  -  l/e$x (x,eH (x) )  .  ( .p) 

Since,  in  general,  $x  =  0(1)  as  e  ->  0  ,  we  have 

Q'(x)  =  0 (1/e )  .  (47') 

The  assumption  in  (12)  was  really  an  anticipation  of  this  condition. 

If,  however,  $x  vanishes  at  some  point,  the  wave  number  in  (47)  is 
undefined  at  that  point.  Our  conclusions  must  be  reconsidered  in  a  neighbor¬ 
hood  of  that  point. 

We  expect  that  there  will  be  a  stagnation  point  on  the  downstream  side  of 
the  body,  presumably  at  the  intersection  of  the  body  and  the  free  surface. 

Such  a  stagnation  point  will  be  located  at  y  =  e/2  ,  as  shown  by  (3) .  We  set 
x  =  x0  at  this  point.* 

At  the  stagnation  point,  we  have  1>X  =  <f>y  =  0  .  If,  separately,  we  require 
that  4>x  =  $  =  0  at  the  stagnation  point,  we  create  precisely  the  condition 

mentioned  above:  The  wave  number  is  undefined.  So  we  consider  more  carefully 
the  behavior  of  <J>  near  (x0,e/2)  .  In  Appendix  A  we  show  that 


e2u0  +  o(e2)  , 

(48a) 

$y(Xg,e/2)  = 

0(e5)  , 

(48b) 

where 

C 

O 

il 

"  C(0)*gxx(Xg,0)  , 

(48c) 

C(0)  = 

body  curvature  at  y  =  0  , 

(48d) 

<(>0  (x,y) 

=  first  term  in  an  e  expansion  of  i(x,y).^ 

(48e) 

*As  in  (I) ,  we 
body  and  the 

shall  also  take 
x  axis.  Since 

x  =  Xg  at  the  downstream  intersection 
the  body  is  assumed  to  be  smooth  and  to 

of  the 

have  a 

vertical  tangent  at  y  =  0  ,  this  practice  should  not  cause  significant  error. 

t4>0 (x,y)  is  the  solution  of  the  "rigid-wall"  or  "double-body"  problem,  in  which 
the  free  surface  is  replaced  by  a  rigid  wall  at  y  =  0  .  See  (I)  for  details. 
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From  these  relationships,  we  now  show  that  the  wave  number  is  G(l/eb)  in  a 
small  region  near  (x0,e/2)  . 

We  substitute  (8)  into  the  free-surface  condition  (5) ,  eliminate  the  term 
that  involve  only  the  naive  expansion,  and  then  determine  the  leading -order 
wavelike  terms.  In  view  of  the  estimates  (48a)  and  (43b) ,  we  can,  to  leading 
order,  drop  all  terms  except  the  following: 

-  —  ?  ~ 

i>  +  <)>  %  0 

y  c-vxvxx 

This  is  directly  comparable  to  the  problem  statement  in  (I) ,  and  it  is  also 
equivalent  to  (29)  above.  However,  it  is  now  being  used  in  a  small  region 
near  the  stagnation  point,  where,  from  (48a),  4>x  =  0(e2)  .  (In  other  words, 
we  now  use  this  condition  everywhere  for  obtaining  the  lowest-order  term  in 
the  wavelike  solution.)  Thus  we  have 

5y  +  0XX’0  (e s)  'V  0  . 

These  two  terms  must  be  of  the  same  order  of  magnitude,  for  otherwise  t> 
would  not  represent  a  wave  motion.  Furthermore,  since  the  potential  satisfies 
the  Laplace  equation,  9/3x  and  d/3y  should  have  similar  order-of -magnitude 
effects  on  5  .  These  requirements  can  be  satisfied  only  if 

k  '  olE'5’  1491 

when  acting  on  4>  ,  which  is  equivalent  to  the  statement  that  wave  number  is 
0(e~5)  .  This  is  valid  only  in  a  region  near  the  stagnation  point. 

From  (48a)  we  can  also  determine  the  order  of  magnitude  of  4>  in  the 
neighborhood  of  the  stagnation  point.  The  complete  potential,  i>  +  $  ,  must 
give  no  normal  velocity  component  on  the  body.  However,  from  (48a)  ,  4>x  = 

0(e2)  at  tne  stagnation  point,  and  so  95/9n  =  0(e2)  there.  In  order  to 
cancel  this,  we  must  have  3>t/3n  =  0(e2)  too  at  the  stagnation  point.  But, 
in  view  of  (49) ,  this  is  possible  only  if 

1  =  C(t7)  (50) 

near  the  stagnation  point.  (Note  that  this  still  does  not  give  us  u  in  (10) 
and  (11).  Those  expressions  are  not  valid  near  the  stagnation  point.) 


The  next  problem  is  to  formulate  and  solve  a  precise  near-field  problem 
for  matching  with  the  outer  solution  from  Section  II. 

We  define  the  near  field  as  a  region  in  which 


x  -  xQ  =  0(e°)  , 

y  -  eH(x)  =  0(e5)  . 


It  is  sometimes  convenient  to  define  new  near-field  variables: 


X  =  (x  -  x  )/eJ 
o 


Y  =  [y  -  eH (x) ]/tJ  j 

We  note  the  following  rules  for  differentiation: 

3  _  - 5  3  -4- 1 ,  3 

3x  '  1  3X  '  M'  (X)3Y  ' 

X  -  e-sX  (53 

3y  e  3Y  " 

Since  we  shall  solve  the  near-field  problem  just  to  one  order  of  magnitude,  we 
shall  have  to  use  only  the  first  term  on  the  right-hand  side  of  the  first 
formula  of  (53) . 

In  the  near  field,  we  expand  the  potential: 


$ (x,y)  =  $ (x,y; e)  +  e7$1(X,Y;E)  +  ... 


From  this  expansion  and  (3),  we  find  that  we  can  also  write: 


H (x)  =  H(x;e)  +  eHHj(X;c)  +  ...  . 


The  free-surface  condition,  (5),  can  now  be  expanded  as  follows: 

0  =  <ty  +  e  ^  { l  ^x  "*■  t  ^xx  t  j  +  •  • .  ]  +  •••} 


—  $>y  +  t^x^xx  +  •  •  •  +  c 


J*lv  +  +  ES,E'3^lvv  +  ' 


to  be  satisfied  on  y  =  eH  =  eH  +  e 5H  j  +  . . .  ,  which  is  equivalent  to  Y  =  Hj  +  . 
In  the  near  field,  $  and  its  derivatives  can  be  evaluated  at  x  =  xQ  ,  y  = 
with  negligible  (higher-order)  error.  Thus,  for  example, 


(5 


4v(x,eH)  =  4x(xQ,e/2)  +  (x-xn)  4VV,  (xn  ,  c/2)  +  e  (H  -  — )  4WV  (xn  ,  c/2)  + 


0 '  xx  1 


2'  yy  o1 


=  <{>x(x0,e/2)  +  ebX$xx  +  U^Hj  +  ebXH  (xQ)  +  ...]  4>yy 


5vd!  -1-  I  5U  4.  c  6  VU  1 

-V  4X  (xQ , e/2)  %  e2u0 

(See  (48a)).  The  free-surface  condition  becomes,  to  leading  order, 


*iy  *  "i/ixx  '  0 


on  Y  =  Hj  +  . . .  . 


The  body  boundary  condition  is  as  follows: 


34_ 

3n 


,34-  ■ 


34  7 L 

r —  +  e  — —  +  . .  .  -  3  . 

3n  3n 


From  the  naive  expansion,  worked  out  in  (I),  we  have 
34 


3n 

34 


=  0 


on  body,  y  <  0  , 


3x  =  “  2eyC(O)40xx(x0,O)  =  2eyu0  on  body,  0  <  y  <  e/2 


Thus  we  require  that 


y34  j 

e7^-  =  0 

3n 


on  body,  y  <  0  , 


?34i 

e /— —  =  -  2eyun  on  body,  0  <  y  <  e/2 

dx  u 


In  terms  of  near-field  variables,  the  last  condition  can  be  written: 


34 1 


3X 


uQ(l+2e4Y)  for  X  =  0  ,  -l/2e4  <  Y  <  0 


<5i 


(59, 


(591 


(6i 


It  would  be  consistent  at  this  point  to  simplify  this  condition  to  S^/SX  = 
-  Uq  and  to  apply  it  in  -<*>  <  Y  <  0  .  However,  this  would  lead  us  to  some 
undefined  integrals,  and  so  we  use  (60)  as  stated  above.  Consistency  can  be 
achieved  later.  We  supplement  (60)  with  the  further  condition: 


34 1 


3x 


=  0 


for  X  =  0  ,  Y  <  -l/2e 


(60 


This  is  consistent  with  (59a)  in  a  small  region  in  which  y  =  0(e)  . 


Finally,  it  should  be  noted  that,  to  leading  order,  4j  satisfies  the 
Laplace  equation: 
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$ i xx  +  ^lyy  =  0  in  the  fluid  r,-,yi'-',u  (6i i 

The  problem  just  formulated  for  4> ;  would  be  straightforward  to  solve 
except  for  one  difficulty:  The  f r<_e-sur face  condition,  (58),  is  to  be  satis¬ 
fied  on  Y  =  Hj  +  ...  .  Usually  in  such  problems  it  is  easy  to  show  that  t.i.c 
boundary  condition  can  be  transferred  to  the  undisturbed  surface  with  negli¬ 
gible  error,  but  such  an  operation  is  not  trivial  in  the  present  i  roblem. 

In  terms  of  near-field  variables,  we  must  note  that  i/lY  -  0(1)  and  H ,  = 
0(1)  ,  and  so  a  simple  Taylor  expansion  cannot  Le  used.  In  terms  of  the 
original  physical  variables,  we  have  itj/ay  -  t>( Jj/c^)  and  the  boundary 

“  f  %  C 

has  to  be  moved  a  distance  eH  -  CH  -  c°Hj  +  ...  =  0{cJ)  .  This  shows  again 

that  the  transfer  is  not  trivial. 


Nevertheless,  it  can  be  carried  out.  Our  demonstration  is  not  rigorous, 
but  it  is  convincing.  We  suppose  first  that  (58)  can  be  applied  on  Y  =  0  . 
Then  we  show  that  the  solution  so  obtained  satisfies  the  same  condition 
applied  on  Y  =  H ^  (to  an  acceptable  accuracy) .  We  then  presume  that  the 
reverse  is  true,  that  is,  that  a  solution  satisfying  (58)  (as  stated)  also 
satisfies  (58)  approximately  when  it  is  imposed  on  Y  =  0  . 


If  (58)  is  satisfied  on  Y  =  0  ,  we  expect  the  solution  to  take  the 

form 


^(X.Y)  =  Fq  (x,y)  +  Fj  (x,y)  exp  {is  (x,y)/e5} 


(62) 


where  the  functions  FQ(x,y)  ,  F^x.y)  ,  and  S(x,y)  all  vary  "slowly"  with 
x  and  y  ,  that  is,  3/3x  and  3/3y  =  0(1)  when  operating  on  these  func¬ 
tions.  (Such  a  result  is  well-known  for  the  corresponding  "wavemaker  prob¬ 
lem."  See  Appendix  R.)  In  fact,  in  the  near  field,  F^x.y)  is  a  constant, 
and 

exp  { is  (x  ,y) /e  =  exp  { iv  (X  +  iY)  }  , 


where 

V  =  1/Ug  , 

and  so  we  see  explicitly  that  S(x,y)  =  v{(x-xQ)  +  i  ly-eii  (x)  ]  }  ,  which  indeed 


varies  slowly  with  x  and  y  .  The  function  Fg(x,y)  represents  a  nonwave¬ 
like  motion,  a  local  effect,  that  decays  with  distance  from  the  body. 

There  is  negligible  difference  whether  we  evaluate  Fg(x,y)  ,  F,(x,y)  , 

and  S(x,y)  on  Y  =  0  or  on  Y  =  ,  and  so  we  need  consider  only  the 

exponential  factor  in  (62) ,  which  does  change  rapidly  with  x  and  y  .  We 
observe  that 

exp  {iS(x,eH)/c5}  =  exp  {i  [3  (x,f:fi)  *-  t:  (il-H)  Sy  (x,eH)  +  .  .  .  ]/c  ‘J  } 

=  exp  {is  (x,EH)/e5}  •  exp  '  iH jSy  (x, ell)  }  •  { 1  +  o  (1)  } 

Thus,  if  we  evaluate  the  exponential  function  on  Y  =  Hg+. . .  (which  is  equi¬ 
valent  to  y  =  eH  )  instead  of  on  Y  =  0  (equivalent  to  y  =  ell  )  ,  we  effec¬ 
tively  multiply  by  a  factor  exp  { iH 1 Sy (x, eH) } 

Now  we  assume  explicitly  that  (62)  satisfies  (58)  on  Y  =  0  ,  that  is. 


(*lv  +  U0$lxy) 


_  (F0v  +  U0F0yy) 


+  U0^2)  tFl  exP  tiS/e5}] 


=  0  . 


The  first  term  on  the  right-hand  side  varies  slowly  in  x  ,  whereas  the  second 
term  varies  very  rapidly.  Then  the  sum  can  equal  zero  only  if  each  term  sepa¬ 
rately  equals  zero.  So  we  have 

F0 Y  +  U0F0XX  =  ° 

s  on  Y  =  0  (  y  =  eH  )  . 

fe  +  uo^  [Fi  exi’ {is/c5}1  =  °  • 

On  the  other  hand,  let  us  evaluate  the  left-hand  side  of  (58)  on  Y  =  Hj+...  : 


(*1Y  +  U0$lxx) 


‘Y^H^. 


=  (Foy  +  uoFoXx^> 


+  exp  {iHjSy  (x,eH)  }  (—  +  Uq— 2}  [Fx  exp  (iS/e5)  ]  I  +  ...  . 

d  9X  'Y=0 

Each  term  on  the  right-hand  side  is  separately  equal  to  zero  (to  the  order  of 
magnitude  considered  here) ,  and  so  we  see  that  (58)  is  approximately  valid  if 
applied  on  Y  =  Hj+...  .  Thus  we  have  shown  the  stated  result:  If  (58)  is 
satisfied  on  Y  =  0  it  is  also  satisfied  approximately  on  Y  =  HJ+...  .  Now 
we  assume  that  the  converse  is  true,  which  means  that  we  can  simplify  our 
boundary-value  problem  by  imposing  (58)  on  Y  =  0  . 


The  solution  will  be  constructed  with  the  followina  G: i en  function: 


G(X,Y;C, n) 


1 


2tt 


log  [  (x-g)  +  (Y-n)']‘  z  +  -—log  +  u  ■.) 


!/2 


+  -  +- 
7T 


dk  k(Y+n)  .  ,v  . .  ..  .-(Yu,) 

e  cos  k(X-  ,)  -  X;  sin  .■  t  X  -  -  )  .  t) 


k  -  v 


M  is  a  constant  to  be  determined.  The  integral  ;s  to  be  interpreted  in  a 
principal-value  sense  (denoted  by  the  bar  through  the  integral  sign).  This 
Green  function  gives  the  potential  at  (>  •  corresponding  to  a  unit  source 

located  at  (£,n)  .  However,  an  audit.! ...  ul  wave  disturbance  has  been  intro¬ 
duced  with  the  M  term.  The  above  Green  function  satisfies  the  Laplace 
equation  and  the  free-surface  condition,  (LG).  It  represents  trie  following 
wave  motion  very  far  away: 


Gx(X,Y;C,n)  ^ 


V  (M 


1) 


,  (Y+n) 
e 


cos  v  (X-i,) 


as  X-c,  ►  ±,v 


(64) 


Using  this  fundamental  solution,  we  express  the  solution  of  our  problem  in 
the  following  form: 


(X,Y) 


0 

vY 

M  (n)  G  (X  ,  Y ;  0 ,  n )  dn  +  De  cos  jX  , 


(65) 


where  p(n)  is  an  unknown  function  to  be  determined  so  that  the  body  boundary 
condition,  (60)  and  (60'),  is  satisfied.  The  extra  term  on  the  right-hand  side 
of  (65)  is  easily  recognized  as  a  solution  of  the  homogeneous  problem,  since  it 
yields  no  contribution  to  ^  at  X  =  0  . 


We  substitute  (c5)  into  the  body  boundary  condition,  obtaining: 


i)$  i  1  vY 

lit  -  2  ^(Y)  --'Me 


/  1  Yn , 
u(n)e  dn 


-  ug(l+2r.4Y)  ,  -1/2:  4  Y  ■:  0  , 


0  , 


(66) 


This  is  a  very  simple  integral  equation  for  , (Y)  .  Tin  solution  can  ohviou  ly 

be  expressed  in  the  form: 


Vi  (Y)  =  M0  (Y)  +  He”''1 

where 


PQ(Y) 


'  - 2u0  (1 +  2.  "Y) 

0  ' 


-!/_•  • 


V  -  1 


(u7 ) 


(On ) 


Substituting  (67)  and  (68)  back  into  (66),  we  obtain  t : -  v  i ;  :■ 
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E  =  - 


4ugM 


1  -  M 


i  - 


The  constants  M  and  D  are  not  yet  known,  but  otherwise  we  know  everything 
in  (65),  and  so  we  have  obtained  a  solution  that  satisfies  (i)  the  Laplace 
equation,  (ii)  the  free-surface  condition,  and  (iii)  the  body  condition. 


Earlier,  we  assumed  that  there  was  a  stagnation  point  at  (Xa,c/2) 
(or  X  =  0  ,  V  =  0  )  .  We  now  find  that,  this  condition  is  sufficient  for 
determining  the  constants  M  and  D  .  1 rem  (54)  and  (93),  we  have 

lb 


-ail 


3y  3y 

From  the  formulation  of  the 


3Y 


(70) 


problem,  we  know  that 


=  -  c  ‘I>2 

9y  X  XX 


’C2  (0)  4>q3  (x0 , 0)  +  o(e5)  , 


the  last  estimate  being  valid  at  the  stagnation  point  (See  Appendix  A) 
must  have 

3,,  2  j 


( 71 ) 
So  we 


d$1  =  +  e3C2(O)^03xx(x0,O) 


(72) 


3  Y  '  ,Jxx  '  u  ■  •  «  '  “xx 

at  X  =  0  ,  Y  =  0  in  order  that  3'l>/3y  ,  as  expressed  in  (70)  ,  may  vanish  at 
the  stagnation  point.  However,  from  the  solution  as  given  in  (65),  we  find 


that  30./3Y  is  undefined  (infinite)  at  this  point  unless  p(0) 


0 


To  see 


this  explicitly,  we  differentiate  (65) : 


34>  i 


3Y 


0 


X=0 


,  ,  3G  (0 ,  Y  ;  0 ,  r| )  ,  .  _  VY 

M  (n)  - ^ — ■LJ—  dn  +  v  D  e 


(73) 


From  the  definition  of  the  Green  function,  (63),  we  have 

k (Y+n) 


3G 

3Y 


X=C=0 


2-n 


y  -  n  Y+n 


n 


dk  k  e 


k  -  v 


^  dk  k  e^1"1 


k  -  v 


as 


0 

Y  +  0 


0 


The  last  expression  behaves  as  -1/uri  as  n  +  0  ,  which  shows  that  the  integral 
in  (73)  diverges  as  Y  ->  0  ,  unless  it  happens  that  p(0)  =  0  .  So  we  now 
impose  such  a  condition  on  p  ,  and  this  then  determines  M  : 


0  =  p  (0)  =  Pq  (0)  +  E 


,  4uqM 

U0  1  -  M 


1  _ 


where  we  have  used  (67),  (68),  and  (69).  Solving  for  M  ,  we  obtain: 


It  is  useful  also  to  note  the  following  consequences; 


E  = 


.0(0) 


U(Y) 


f  2u0{ev  -  1  - 
VY 


v.  “UUe 


Y  0 
y  <  -l /2c1' 


r/'  > 


The  above  choice  of  V  guarantee..  t».at  '  , 
nation  point,  but  it  does  not  yet  lead  to  (72).  In 
is  satisfied,  we  must  evaluate  the  integral  in  (72) 
the  value  of  D  appropr iately .  The  evaluation  of 
in  Appendix  <:.  From  these  results,  we  obtain: 

D  =  c3u0*0xx<x0'°>/  =  :  U0'0xx(x0'J>  = 


remains  finite  at  the 
order  to  ensure  that 
as  Y  >  0  and  then 
the  integrals  will  be 


<d)<iCxx(x0,t!) 


s  tag  - 
(72) 
chooss 
;  ound 


(77) 


In  summary,  we  note  that  the  solution  in  the  inner  region  is  expanded  as 
in  (54),  <j>  being  the  naive  solution  that  was  worked  out  in  (I).  The  next 

term,  ,  is  obtained  in  the  form  given  in  (65),  with  t.  (Y)  given  in  (76) 

and  D  in  (77).  This  completes  the  inner  solution  for  our  [urposes. 


Before  matching  this  inner  solution  with  the  outer  solution,  it  is  worth¬ 
while  to  comment  on  some  unusual  aspects  of  the  foregoing  analysis: 


(i)  Normally  in  using  matched  asymptotic  expansion  ,  one  co tains  inner 
and  outer  expansions  each  of  which  is  nonunique  in  some  way.  Matching  of  the 
two  expansions  then  removes  the  nonuniqueness  in  each .  However,  in  our  prob¬ 
lem,  we  have  obtained  a  complete  (unique)  near-field  solution,  at  least  to 
leading  order  of  magnitude.  This  appears  at  first  sight  to  eliminate  the 
possibility  of  satisfying  a  radiation  condition  in  the  far  field.  In  fact, 
this  is  precisely  what  happens,  and  fortunately  too,  for  there  is  no  radiation 
condition  possible  in  the  far  field.  In  steady -motion  problems,  we  usually 
specify  that  there  should  be  no  waves  upstream;  this  is  an  adequate  radiation 
condition  in  two  dimensions.  In  our  problem,  however,  the  downstream  waves 
appear  only  in  a  surface  layer  of  vanishing  thickness,  and  these  waves  are 
effectively  isolated  from  the  upstream  fluid  region.  One  might  say  that  there 


is  no  upstream  region  at  all,  at  least  insofar  as  it  might  affect  the  waves 


downstream  of  the  body.  Then  one  must  seek  an  alternative  condition  to  make 
the  solution  unique.  As  described  above,  we  have  used  the  solution  behavior 
at  the  stagnation  point  to  provide  such  a  condition. 

(ii)  The  Green  function  introduced  in  (63)  is  rather  unusual  because  of 
the  presence  of  the  last  term.  The  contribution  of  this  term  to  1 .  is 
interesting.  Effectively,  it  produces  a  second  solution  of  the  homogeneous 
problem.  The  one  obvious  solution  of  the  homogeneous  problem  is  the  last 
term  in  (65).  But  one  can  construct  ethers  if  some  degree  of  singul  irity  is 
tolerated  at  the  origin.  The  M  term  in  the  Green  function  produces  a  solu¬ 
tion  of  the  homogeneous  problem  with  precisely  the  singular  character  at  the 
origin  that  is  needed  to  cancel  the  singular  velocity  that  would  otherwise 
arise  from  the  source  distribution.  This  is  what  was  accomplished  in  deriving 
the  value  of  M  in  (74) . 


IV.  MATCHING  OF  INNER  ANil  OUTER  SOLUTIONS 


We  now  use  the  inner  solution,  which  is  completely  known,  to  Jet  :rmi  n-.- 
the  unknown  amplitude  and  phase  of  the  waves  ii.  the  outer  region.  From  (8) , 
(10),  and  (38),  the  outer  solution  can  be  written: 

4>(x,y)  =  $(x,y;e)  +  ea+1  exp  [  (Y  -  H) /1>X]  Re  {a (x,y)  exp  [if< (x)  j  }  .  t"7 

The  rapidly  varying  phase  function,  0(x)  ,  was  given  in  (37);  the  complex 

amplitude  function,  A(x,y)  ,  was  partially  determined: 


A(x,0)  =  AQ/4>X  (x,  uH  (x)  )  , 


(7  ) 


where  A^  is  a  complex  constant.  All  that  remains  to  be  determined  is  A 
and  a  .  We  accomplish  this  by  matching  4>x(x,eH)  in  the  two  regions. 

For  the  moment,  let  us  omit  the  Re  notation  in  (78),  just  as  we  did 
earlier.  We  write  out  <Jx(x,y)  : 


4x(x,y)  =  4x(x,y,-e)  +  e 


_a+l 


,  H’  2(Y-H)'4xx 

Ax(x,y)  +A(x,y){-^ - -3  -  ■  +  i0‘  (x)  j  j 

*5  ^x  ’ 

•  exp  [  (Y-H)/$x]  •  exp  [i<3  (x)  ]  . 


On  y  =  cH  (or  Y  =  H  *),  this  simplifies: 

4>x(x,eH)  =  4>x(x,eH;e)  +  ea+1  |"ax(x,0)  +  A(x,0)  {-iy+ i0’ (x)  } 


+  ... 


(80) 

exp  [iO(x)J 
(8C) 


-4 

» 


n 


Now  we  change  over  to  near-field  variables,  as  defined  in  (52) .  We  note  that 

1 


O’  (x)  = - =m— ^ — 

e4>x  (x, £H) 


'V  - 


^ — y  as  X  -*■  0  ; 


e  "u 


0  (x)  =  -  - 

£ 


dC 


e4u^[l +  0(e3) ] 


x~  x0 

c  5. .2 

e  Ug 


*Note  the  difference  in  the  definition  of  Y  in  the  inner  and  outer  regions: 

Y  =  y/c  in  the  outer  region  (see  (13)),  whereas  Y=  [y-£H(x))/e5  in  the  inner 
region  (see  (52)). 


A 

■M 
■  1 
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The  estimate  (48a)  has  been  used  here.  Now  we  keep  just  the  leading  term  in 
(80'),  with  X  fixed: 


<J>x(x,eH)  n,  $x(x,eH;e)  -  i (A0/uq)  ea-6  e  lX//u  0 
This  is  the  result  that  we  will  match  to  the  inner  solution. 


(81) 


The  inner  solution  was  given  in  (54) ,  (65) ,  (76) ,  and  (77) .  Let  us  find 
the  derivative  that  matches  (81): 

<i>x(x,y)  =  $x(x,y,-e)  +  e7[^1x^x  +  ^iyYx)  +  ... 

=  4x(x,y;e)  +  c2*lx  +  ...  - 

We  must  evaluate  this  as  X  ->  °°  in  order  to  perform  the  matching.  Introducing 
(64)  into  (65) ,  differentiating  with  respect  to  X  ,  and  evaluating  the  result¬ 
ing  integral,  we  find  that 

4>lx  n,  -v{4el+u0  cos  vX  +  Dsinvx}  eVY  t  -  e5c2  (0) (Jjq3^ (xQ , 0)  sin  X/u? 

+  ...  .  (82) 


In  (82),  the  term  containing  cos  vX  comes  from  the  integral  in  (65),  and  the 
term  containing  sin  vX  comes  from  the  homogeneous  solution  in  (65) .  As  is 
apparent  here,  the  latter  is  of  lower  order  of  magnitude  and  thus  dominates  in 
(82)  . 


Now  we  match  these  near-field  results  with  the  real  part  of  (81) ,  the 
result  being  that 

AQ  =  -  C5(O)^h(x0,O)  .  (83) 

We  also  find  that 

a  =  11  .  (84) 


We  can  substitute  back  into  the  formulas  for  4>  or  for  4>x 
The  latter,  in  particular,  gives  us: 

, ,  (0)  <j> q*xx  (Xq  ,  °) 

$x(x,eH(x))  =  3>x(x,eH;e)  +  xe11 - =2  - exp  - 

^x 

This  is  our  final  form  of  the  solution  in  terms  of  4>  . 


in  the  far  field. 


*  \ 
if  ♦ . 
$1  I 


(85) 


V. 


WAVE  RESISTANCE 


Exact  nonlinear  formulas  are  readily  available  for  computing  wave  resis¬ 
tance  if  the  fluid  velocity  and  free-surface  disturbance  are  known .  From 
Wehausen  and  Laitone  (i960),  Equation  (8.6),  for  example,  we  have: 

^eH(x) 

R  =  ~  P  U“L  [  dy  {^y(x,y)  -  (4>x  (x,y)  -  1]  2|  t  jpgL^¥(x)  ,  (86) 

^  J  * 

— oo 

where  R  is  the  wave  resistance,  and  all  other  quantities  are  as  previously 
defined.  The  right-hand  side  can  be  computed  at  any  x  downstream  of  the 
body,  but  we  simplify  the  task  by  letting  x  -*■  “  . 

The  general  far-field  expression  for  4>x(x,y)  has  been  given  in  (80).  As 
x  -*■  +oo  ,  the  following  approximations  can  all  be  used: 

$x  -  1  #  ^xx  ~  0  ,  $y  =  0  ; 

A (x,y )  -  A0  ,  Ax  =  0  ; 

(87) 

H  -  0  ,  H1  =  0  ; 

O'  (x)  =  -  1/e  . 

Noting  that  AQ  is  a  real  constant  and  taking  the  appropriate  real  parts  in 
(80)  (see  (78)),  we  obtain  the  asymptotic  estimate, 

~  1  +  e 1 1  AQey^E  sin  [0(x)]  ,  (88a) 

valid  as  x  +»  .  Similarly  we  find  that 

^y(x,y)  -  e^Ajje^^6  cos  [0(x)]  .  (88b) 

Moreover,  from  (9),  (11),  (39),  and  (40),  we  also  have: 

H(x)  =  -  e^AgSin  [0(x)]  .  (88c) 


Carrying  out  the  integration  in  (86) ,  with  the  upper  limit  of  the  integral 
replaced  (consistently)  by  zero,  we  obtain  finally  that 
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where,  as  before,  C(0)  is  the  body  curvature  at  the  ■  ntersecti of  the  body 
and  the  undisturbed  free  surface  and  ‘Jq(x,v)  is  the  first  term  in  the  naive 
expansion  (the  potential  for  the  double-body  flow) .  Cw  is  the  wave -resistant' 
coefficient  for  the  2-D  body. 

It  can  hardly  be  surprising  that  the  wave  resistance  depends  op  the  shape 
of  the  body  only  near  the  free  surface,  since  the  wave  motion  occurs  entirely 
in  a  very  thin  layer,  in  fact,  in  a  vanishingly  thin  layer,  near  the  level  of 
the  free  surface,  and  the  wave  motion  i  induced  by  the  peculiar  nature*  of  t  ’:><• 
streaming  flow  on  the  body  near  this  level  (see  (I)).  Furthermore,  we  assumed 
that  the  body  has  a  vertical  tangent  here,  and  so  it  is  also  reasonable  to 
expect  that  the  curvature  should  have  a  dominant  effect  in  creating  waves. 

If  the  curvature  were  .*,.-.ro  at  y  =  0  ,  we  would  expect  the  generation  of 
waves  to  depend  primarily  on  higher-order  derivatives  of  the  body  shape  in 
this  same  region.  Presumably  wave  amplitude  and  wave  resistance  would  then  be 
of  even  higher  order  than  in  the  case  presented. 

We  can  only  speculate  now  what  would  be  the  result  if  the  body  contour 
were  analytically  straight  in  some  finite  neighborhood  of  y  =  0  .  Our  specu¬ 
lation  is  that  waves  and  wave  resistance  would  be  small  of  exponential  order 
as  in  the  case  of  a  submerged  body . 
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APPENDIX  A. 


STAGNATION- PC  I  NT  CONDITIONS  ON  4> 


Here  we  obtain  the  estimate  given  in  (48 c)  and  (71) , 


From  (I),  from  the  definition  of  the  naive  expansion, 


*x  %  *0X  +  ^lx  +  •••  • 

On  the  body  very  close  to  the  free  surface. 


If  the  body  shape  has  continuous  curvature  and  a  vertical  tangent  at  y  -  <'  , 
we  have,  from  Equation  (42)  of  (I), 

*ox  =  o(e3)  . 

From  Equation  (62)  of  (I), 

<hv|  -  ~  2  y  C(O)<f>0  (x0,O)  for  0  <  y  <  e/2  , 

1  body 

where  C  (0 )  and  <J>Q(x,y)  are  given  by  (48d)  and  (48e)  .  Putting  these  results 
together  gives  us: 

5x(xq,§0  =  -  e2C (0) 4>0xx (xQ ,0)  ,  (A.l) 

the  relationship  given  in  (48a) . 


To  obtain  (71)  (and  also  (48b)),  we  start  from  the  free-surface  condition 
(5)  ,  which  is  still  valid  (by  definition)  if  $  is  replaced  by  4>  .  Thus 


+  2$  $  +  $^5  } 

xx  x  y  xy  y  yy1 


on  y  =  eH(x) 


On  the  right-hand  side,  the  first  term  is  the  lowest-order  term,  and  so,  sub¬ 
stituting  from  (A.l),  we  obtain: 

Vxo'f}  =  "  e5c2(O)<j>03xx(x0,o)  ' 

the  desired  result,  (71). 


APPENDIX  B.  COMPONENTS  OF  THE  INNER  SOLUTION 


In  (62),  the  inner  solution  was  written  as  the  sum  of  two  parts: 

>t’,(X,Y)  =  Fq  (x,y)  t  Fj  (x,y)  exp  {is  (x,y)/t  f'}  ,  (62) 

where  it  was  stated  that  F-  (x,y)  ,  F] (x,y)  ,  and  S(x,y)  all  vary  "slowly" 

with  x  and  y  ,  although  this  expression  represents  waves  of  very  small 
wavelength. 

Such  statements  can  be  verified  directly  in  the  analogous  "wave-maker 
problem: " 


$xx  +  *yy  =  0  in  |  *  <  o  i 


(6.1) 


<C x  =  f  (y )  on  x  =  0  ,  y  <  0  , 


(B.2) 


vc}>  —  <Jiy  =  0  on  y  =  0,x>0 


(B.3) 


This  is  the  mathematical  statement  of  the  problem  if  the  vertical  wall,  x  =  0  , 
moves  with  normal  velocity  component  f(y)  exp  (icot)  ,  where  w  is  the  radian 
frequency  of  the  motion  and  v  =  w2/g  is  the  corresponding  wave  number.  The 
potential  for  the  problem  is  $(x,y)  exp  (iiot)  .  Of  course,  a  radiation  condi¬ 
tion  must  also  be  imposed.  This  problem  is  very  similar  to  the  near-field 
problem  considered  in  Section  III. 

The  wavemaker  problem  has  been  studied  by  many  people,  and  we  can  write 
down  the  solution  directly: 


(x,y)  =  $  (x,y)  +  $  (x,y)  , 


(H.4) 


where 


<{>(x,y)  =  e  [Ar,  sin  vx  +  Bcosvx)  , 


(1J.5) 


(x,y)  =  dk  A(k)c  (k  cos  ky  +  v  sin  ky )  . 


(B  .6) 


The  term  in  (B.5)  containing  B  is  the  solution  of  the  homogeneous  problem, 
and  so  B  is  arbitrary.  The  constant  An  in  (B.5)  and  the  function  A(k) 
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in  (B.6)  are  given  by  the  following: 


A,,  =  2  I  dy  f  (y)  e 


(B.7) 


A(k)  =  -  — — 5 - dyf(y)  (k  cos  ky  +  v  sin  ky) 

Tik(k2+v2)  ! 


( r. .  8 ) 


In  the  above  problem,  we  allow  v  (or  o>  )  to  become  very  large,  it 
is  immediately  obvious  that  <j>x  and  <j>y  are  0(vf>)  .  We  can  express  this 
symbolically  by  writing 


—  ,  7-  =  0 (v)  . 

ux  3y 


( }  ■ .  9 ) 


Since  if>  represents  waves  of  wavenumber  v  ,  this  is  all  rather  obvious. 


What  is  not  so  obvious  is  that 


3x  '  cy 


=  0(1) 


(B. 10) 


when  they  operate  on  <f>  (x,y)  ,  even  if  v  ->■  00  .  To  show  this,  substitute  (B.8) 
into  (s.6),  interchange  order  of  integration,  and  write  the  potential  in  terms 
of  a  function  of  a  complex  variable  z  =  x + iy  : 


(x,y )  =  -  —  Re  dn  f  (n) 

7T 


00  , 

f  dk  e 

J  k  (k-i\ 


(k  cos  kn  +  vsinkn) 


=  -  Re 

IT 


J  dn  f  (n)  (  log  -  2e 


•iv  (z+in) 


Ej (-iv (z+in) ))  , 


where 


clt  G 

E1  (u)  =  - - -  ,  the  exponential  integral 


As  v  -*■  “  ,  note  that 


E,  (vu)  'v  - 

1  vu 


for  |arg  u]  <  3tt/2 


Using  this  asymptotic  estimate,  we  now  find  that,  as  v  <*>  , 


♦*-i*y  =  V  dnf(n)  bfe  -  0(1/v)  ; 
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^xx  ”  1<^xy 


dn  f  (n) 


(z-ir\)2  (z+tn) 


—  +  O  (1/ vi)  ; 


etc.  This  is  the  result  indicated  by  (R.10).  It  means  that  f(x,y)  varies 
"slowly"  even  if  the  wavenumber  becomes  asymptotically  large. 

In  the  problem  in  Section  III,  we  had  the  free- surf ace  condition 


V(J)y  +  <f>xx  =  0  on  y  =  0  , 


(B.ll) 


instead  of  (8-3).  The  solution  can  be  written  as  in  (B.4),  with  4>(x,y)  giver 
by  (B.5)  and 


(x,y)  =  dkA(k)  e  (v  cosky  -  k  sin  ky)  . 


In  this  problem,  one  can  show  that 


A0  =  2  dy  f  (y)  {l  -  eVy }  ? 


(!' .  12 ) 


(8.13) 


A  (k)  = 


irk  (k2+v 


f° 

dy  f  (y)  { v  (1  -  cos  ky)  +  k  sin  ky}  .  (8 . 14) 


Then  it  is  straightforward  to  show  results  identical  to  those  obtained  above 
for  tne  wavemaker  problem,  as  summarized  in  (8.9)  and  (8.10). 


APPENDIX  C. 


EVALUATION  OF  SEVERAL  INTEGRALS 


We  chose  the  arbitrary  constant  M  in  the  Green  function  (see  (63))  to 
ensure  that  9$^/9Y  would  be  bounded  at  X  =  0  even  as  Y  -*■  0  .  The  result 
was  given  in  (74).  This  left  the  constant  D  ,  which  first  appeared  in  (65), 
as  a  still  undetermined  quantity.  Its  value  was  to  be  found  by  substituting 
from  (73)  into  (72)  .  To  carry  this  out,  we  have  to  evaluate  the  integral  in 
(73)  for  X  =  Y  =  0  .  Its  value  turns  out  to  be  negligible,  but  this  result 
is  not  obvious,  and  so  we  derive  it  here. 


From  (73)  and  (72),  witn  Y  =  0  in  the  former,  we  have 


e  Voxx(V0)  -  vD  =  I  dnu(?)  37 


(r  .1) 


X=Y=0 

U=o 


V  =  1/Ug  . 

,  )  Following  (73) , 

3G 

.  r00,,  .  kn 

1  dk  k  e 

3Y  X=Y=0 

TT  T  k  -  V 

0 

5=0 

p(n)  was 

given  in  (76).  We 

Q 

> 

1 

O 

o 

X 

o  fO  f 

2u0  vn 

=  — dn  e 

71  i. 

kn 


T  J  J,*-v 

—co  0 

2u0  f°  4 

— -  j  dn  (1+  2e4n) 

7T  J 

-1/2e4 


dk  k  e 
k  -  v 


4uq£^ 


7TV  I 

< 

4une4 


(C.2) 


{ Y  +  log  — — r-  -  — +  . . . }  as  v  ■+  oo  ,  (C .  3 ) 


2e  v 


where  y  is  Euler's  constant.  The  expression  in  (C.2)  is  found  after  some 
manipulation  of  the  preceding  line,  and  the  estimate  in  (C.3)  follows  from 
standard  asymptotic  estimates  of  exponential  integrals. 

The  right-hand  side  of  (C.3)  is  0(e4  log  4)  ,  whereas  the  first  term  on 
the  left-hand  side  is  0(e3)  .  So  we  conclude  that 


D  =  e  U0^0XX<X0-0)/V  * 


(C.4) 


as  stated  in  (77) . 
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